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Abstract 

We present a classical analysis of the transient response of Josephson junctions perturbed by 
microwaves and thermal fluctuations. The results include a specific low frequency modulation in 
phase and amplitude behavior of a junction in its zero-voltage state. This transient modulation 
frequency is linked directly to an observed variation in the probability for the system to switch 
to its non-zero voltage state. Complementing previous work on linking classical analysis to the 
experimental observations of Rabi-oscillations, this expanded perturbation method also provides 
closed form analytical results for attenuation of the modulations and the Rabi-type oscillation 
frequency. Results of perturbation analysis are compared directly (and quantitatively) to numerical 
simulations of the classical model as well as published experimental data, suggesting that transients 
to phase-locking are closely related to the observed oscillations. 

PACS numbers: 74.50. +r,03.67.Lx,85.25.Cp 
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I. INTRODUCTION 



Much experimental attention has recently been given to the topic of Rabi-oscillations 
in microwave irradiated Josephson systems at very low temperatures 3, 3, U, ^|. A 
common motivation for the studies is the proposed macroscopic quantum behavior outlined 
in Ref. which has opened the way for potential applications of Josephson technology in 
quantum information processing [t|. This is proposed by operating a Josephson junction 
in its zero-voltage state at temperatures below the quantum transition temperature, while 
manipulating the possible energy states with the application of commensurate microwave 
frequencies. The significance of Rabi-oscillations is many fold, and includes i) a direct 
connection to a well known quantum mechanical concept in perturbed atomic systems, and 
it) a method by which the control of quantum states of a system can be evaluated. Thus, 
it is of significant importance to understand the nature of these observations in order for us 
to evaluate how to interpret the Josephson system under investigation. 

Recent classical analyses of Josephson junctions, perturbed by microwaves and low tem- 
perature thermal fluctuations, have revealed that key signatures of microwave- induced multi- 
peaked switching distributions J^, 0, [h], 1]], used to display quantum mechanical features 
of Josephson junctions, have direct classical analogs, which may obscure the interpretation 
of the observations 112 



13l Il4 Il5| . Inspired by the work on slowly modulated transients to 
phase-locking in Josephson systems by Lomdahl and Samuelsen Q], it has most recently 
been demonstrated (l2| that also Rabi-oscillations have a classical analogue in microwave 
perturbed Josephson junctions, providing a system response very similar to reported obser- 
vations under the same conditions. We therefore denote this classical phenomenon Rabi-type 
oscillations. It is the aim of this paper to complement and expand on the analysis provided 



in Ref. [17| in order to extract more detailed information from the classical model as it re- 
lates to observations of Rabi-oscillations. We further demonstrate the validity of the analysis 
through direct quantitative comparisons with the observations of Ref. j^j], in which enough 
information is provided to conduct direct quantitative comparisons with our theory. Specif- 
ically, we present an analysis based on direct perturbations in the dynamical variables. This 
approach results in closed form results that include dissipaion, which is of direct importance 
for experimental measurements of, e.g, coherence time. 

We will start section II by providing the system equations that govern the classical model 
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as it relates to specific reported measurements of Rabi-oscillations. Next, in section III, we 
present a new, more complete, perturbation analysis of the slow transient modulations to 
phase- locking, and provide specific analytical results for (Rabi) oscillation frequency and 
attenuation as a function of the system parameters. The theory is numerically validated in 
section IV. Direct observations of Rabi-type oscillations in the classical model are produced 
in section V through numerical simulations of a system parameterized to mimic reported 
experiments, and close agreement is found between these simulations, the developed theory, 
and the reported experiments. We conclude the paper in section VI. 

II. CLASSICAL MODEL 

A normalized classical equation for a Josephson junction can be written [l^ 



where (p is the difference between the phases of the quantum mechanical wave functions 
defining the junction, rj represents the dc bias current, and e s (t) and u s represent microwave 
current amplitude and frequency, respectively. All currents are normalized to the critical 
current J c , and time is measured in units of the inverse plasma frequency (jJq, where Uq = 
2eI c /hC = 27rJ c /$oC, C being the capacitance of the junction and $o — h/2e the flux 
quantum. The phase of the applied microwave field is 6 S . Tunneling of quasiparticles is 
represented by the dissipative term, where a = hu /2eRI c is given by the shunt resistance 
R, and the accompanying thermal fluctuations are defined by the dissipation-fluctuation 
relationship 19] 



T being the temperature and Hj = 7 c $ /2vr is the characteristic Josephson energy. A 
current pulse for probing the state of the system is represented by e p (t). 

From the above normalizations, we define the normalized energy from the time indepen- 
dent perturbation terms in Equation 



<p + cap + sin ip = 77 + e s (t) sm(u s t + 6 S ) + e p (t) + n(t) , 



(1) 




(2) 
(3) 



H = —ip 2 + 1 — cos ip — Tjip . 



(4) 
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The undamped bias current dependent plasma frequency is given by 



1 -r/ 2 



(5) 



Within this model, we will make analysis and conduct simulations according to experiments 
reported in Ref. pj. System parameters, such as 77, a, lu s , e s (t), and e p (i) are matched as 
closely as possible to reported values, and Rabi-type oscillations are observed statistically 
through simulations of the distribution of probe field induced switching from the zero- voltage 
state as a function of applied microwave field e s (t) sm(uj s t + 6 S ). Microwave and probe 
perturbations have the form sketched in Figure ^ We assume the phase 8 S of the microwave 
field is random for each switching event. Notice the difference in probe field between this 
presentation (along with Ref. |5[) compared to Refs. , where the probe was a 

microwave field with frequency slightly lower than uj s . 
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FIG. 1: Schematic of signaling protocol, which is adopted from Ref. Dotted line represents the 
dc bias rj which is constant throughout each run. The normalized time ineterval Ato represents a 
91 ns duration during which no ac signal is present to allow the junction to achieve a steady-state. 
At s represents the duration of the ac signal. This value is varied so as to observe the temporal 
Rabi-type oscillations. Values for the signal timings are At a = 1.5nsujQ, Atb = 2.5nsu>o, and 
At c = 2.0nso>o. 
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FIG. 2: Comparison between numerical results of Equation JQ) (markers) and the monochromatic 
ansatz Equation (JSJ) (curves) for a = e s (t) = e p (t) = T = 0. Results for 77 = 0.904607 are displayed 
with o and solid curves, while results for rj = 0.98 are displayed with • and dashed curves, (a) 
Anharmonic relationship between amplitude A and frequency uo r Equation (jHJ); (b) amplitude A and 
average phase (fo Equation (J7J) ; and (c) amplitude A and normalized energy H — Hq Equations @ 
and JTUJ). 

III. PERTURBATION ANALYSIS 



Following the analysis of Refs 



.0 



14. 



201 ] . we will use the monochromatic ansatz for a 



solution of the unperturbed Equation ((TJ) (a = e s (t) = T = e p (t) = 0) 

if = ip + ip = (po + A sm(ujt + 9) , 



(6) 



where is the average phase, A is an oscillation amplitude (which is not necessarily small), 
and io is the natural frequency. Inserting this ansatz into Equation we obtain the 
anharmonic relationships between the parameters 



J (A) sin<^ = V 
2Ji(A) 



V 



MA), 



(7) 
(8) 



where J n is the nth order Bessel function of the first kind j21|. Notice that the ansatz 
dictates Jq(A) > rj as a condition for existence of a zero- voltage state. Also notice that 



w r < and that uo r — > uoj for A — > 0. The average system energy of the monochromatic 
ansatz is found by inserting Equations © and (J7J) into (jlj), then average over the period 
given by (JSJ), 

(H)*. = -A 2 u 2 + 1- J (A) cos (po-r)(p , (9) 
where the minimum energy Hq for a given r/ is 



Hq = 1 — y 1 — r] 2 — r^sin 1 r/ . (10) 

The applicability of this very simple monochromatic ansatz © is illustrated in Figure 121 
Here we have shown the average phase ipo and oscillation frequency u r as a function of 
oscillation amplitude A for two different values of dc bias rj. Data are obtained from direct 
numerical simulations of Equation (markers) and from Equations (|7|)-()10|) (curves). The 
agreement between the results of the ansatz and the simulations is very close, except at the 
highest amplitudes (when J (A) — > rj), where higher harmonics of the dynamics contribute 
significant components to the dynamics. We use this ansatz to develop a perturbation theory 
for the system response to the application of microwaves and damping. 

For T = e p (t) = we now separate the parameters of the ansatz into steady state and 
small time dependent deviations; A = A + 5A, 9 = 9 + 59, and ip = ip + 8(p , where 
\8A\ < 1, \59\ < 1, and \5(f \ < 1. Inserting Equation © into (QJ for e p (t) = T = and 
e s (t) = Q(t)e K (9 being Heavyside's step function), yields the steady state phase- locked 



relationships 



el 



A 2 (K 2 - u 2 r ) 2 + a 2 u 2 ) (11) 



auj s 



tan(0 - 9 S ) = , (12) 

and the linearized expressions for the small deviations from steady state: 

S(p + aSip + Jo (A) cos (fo 6<po = Ji(A) sin<^ SA (13) 
SA + a(5A — uj s A59) - 2uj s A59 + (j (A) - J 2 {A)) cos^ - ^ 2 SA = 2J X (A) sin^ 5(p 

(14) 

A59 + a(A59 + u s 5A) + 2u s 5A + [(j (A) + J 2 (A) cos ^ - oo 2 ] A 59 = . (15) 

Equation (jT3j) represents the small amplitude, slowly varying terms in (JTJ), while (fT%Jl and (fT5)l 
represent small amplitude terms oscillating with uj s . Evaluating a solution to Equations (HHJ)- 
(fl3|) is not simple, but making the assumption that all three variables oscillate slowly with 



the same frequency Qr and decay with the same attenuation j3, we can write the simple 
relationships: SA = exp(zf2t), A58 = k,5A, and <5<£>o = l^A, with Q = Qr + ij3. This results 
in the following polynomial 



Q 6 - 3iaQ 5 
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-a 



(14 



-a 



2a A 



n 3 + 
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— a a a,A + ao 
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uj s a + Ti = a uj s Jq(A) cos<^ + r2 



(16) 



where the parameters <2j and Tj are given in Equations (fTH|) - (|2*2*j) below. The polynomial 
can be simplified considerably by inserting Q = Qr + ij3, and realizing that f3 — -^ol is the 
complex part of the roots (for all roots) with a non-zero real part Qr; i.e., for underdamped 
modulations. Notice that dissipative systems (a > 0) will naturally imply that Q 2 R — > + 
for a nonzero threshold value of the microwave amplitude e s > With (3 = \a the 

(real) modulation frequency Q R is determined by the following real polynomial: 



Q 6 R + a 4 Q R + a 2 Q 2 R + a = 



(17) 



where 
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2J (A) cos^o) + 2(J (A) cosy? + u; 2 ) J (A) cosy2 
- 2J (A) cos(^o) Jo(A) cos<^o + 2J 2 (A) sin 2 y? • 



(18) 
(19) 
(20) 

(21) 
(22) 



The steady-state resonance frequency u r is given by Equation © as a function of A. Gen- 
erally speaking, the three roots Q R are real and positive. One is located near (2ui s ) 2 , one 
near u 2 , and one much smaller than u 2 . It is the latter we are interested in when studying 
transients and modulations to phase-locking, and it can be very well approximated by either 
of the following approximations, since the normalized frequency is small: 



A further simplification to the coefficients in the polynomial can be produced by noticing 
that the damping parameter a in relevant Josephson experiments is very small, allowing for 
omission of several terms of high order in a in Equations (|18 |) -(|2(J |) . A specific limit of the 
above solution can be derived for small oscillation (A) and microwave (e s ) amplitudes for 
a = and u>f = 1 — rf 

n R « M (25) 

V 1 1 

16 c; s 3 16 (l-7? 2 )f 

^^^Jli^^^o. (27) 

4 w| 4 (l-r/ 2 )cu s 12 

The resonant choice uj^ « 1 — 7/ 2 is consistent with the experimental studies ||. In fact, 
measurements are usually conducted for cj s very close to the linear resonance frequency (see, 



e.g., Ref. 2) 14)- However, if ^ 1 — r] 2 , then A oc e s , and the relationship between Q 2 R 
and e 2 s becomes linear for small values of A and e s . Further, if uj^ < 1 — rj 2 , A may become 
a multi- valued function of e s , corresponding to several resonant energy states of the system 
for the same parameter values [171 ] . 

The solution to Equation ([17jl. which can be given in closed form as a solution to a third 
order polynomial, is a complementary, more convenient, approach to the one presented in 
Ref. [171 . It provides the important additional information regarding the attenuation (3 of 
the transient modulations and it provides explicit expressions for the modulation frequency 
as a function of all system parameters. The two different perturbation methods agree qual- 
itatively, and are quantitatively similar. The main difference between the two approaches is 
that we have here considered phase 9 and amplitude A as harmonically varying variables, 
instead of using the energy balance approach to transients |l6j |. which treats the total energy 
of the system as the linearized dynamical variable. Notice that the linearized Equations ([13)1 - 
(115)1 provide information about frequency, attenuation, and internal phase-relationships, but 
not about the specific magnitude of the deviations. However, the magnitude can be quan- 
titatively estimated from the following reasoning. For a simulation, in which a Josephson 
junction is described by Equation (JTJ) with e p (t) = T = (or at least very low temperature) 
and e 3 (t) = Q(t)e s , the system is at rest ip = <p = sin -1 77 for t < (A = 0). The onset of 
the microwave field at t — will therefore, within the harmonic approximation, produce an 
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envelope function of oscillation given by 22 



A — A + 8A — { , ' - . (28) 

A(l- e-f***n R t) , t > 

The two other modulated variables, 5(fo = j5A and AS6 = k,5A, can then be found from the 
coefficients 

7 = -Ji(A)sm<p 

f2 2 — iaVt — Jq (A) cos (fo 
auj s + 2iuj s VL 



VI 2 - iaQ + u 2 s - uo 2 ' 

The modulated system energy (for e p (t) = 0) can be calculated from Equation © as 

(H)2n = H + 5H , (31) 

uj s 

where H is the steady state energy of the phase-locked state in Equations (fTTj) and (|T2"|). 
and 5iY is the transient modulation. The average phase can be calculated either as a 
function of A = A + 8A as given by Equations and (|28j) . or it can be extracted from the 
theory through Equation (j2H|l . The two results are very similar, but since the latter is a first 
order approximation in the perturbation theory, it will have a second order error at t = 0. 
Using Equation (JJJ), <£>o(£) wm satisfy the true value both at t — and for £ — » oo. 



IV. SIMULATION OF TRANSIENTS 



Numerical validation of the expressions for the transient modulation can be directly 
acquired from simulating Equation (p^ for e p (t) = T = and e s (t) = e s Q(t) for different 
values of = 1 — rf and a. Simulations are conducted by choosing r\ < 1 and a > 0, 
then initiating the system in the static state (<p(0),(p(0)) = (sin _1 ?7,0). At t — 0, the 
microwave field switches on, and we measure the phase </?(£) i, the average phase ^o(^) = 
(<fi(t))22L, and the system energy H(t) as given by Equation (jlj). Figure El shows a typical 
comparison between the analysis of the previous section and direct numerical simulations. 
We have chosen system parameters a = 0.00151477, rj = 0.904706, e s = 0.00108 ■ O(t), and 
uj s = v^l — rj 2 = 0.6527147 (inspired by Ref. [5j). We can clearly see that the transient in 
phase and amplitude is a very slow modulation, and that this modulation provides a slow 
transient oscillation in the system energy. Notice that the phase-locked frequency is much 
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higher than the depicted modulation, <C 00 s . We have shown results of the theory on 
the left (Figures and b) and the corresponding simulation results on the right (Figures Et 
and d). The quantitative agreement between the results is obvious. We do see a slightly 
larger predicted oscillation amplitude compared to the simulated. We also notice a slightly 
higher predicted modulation frequency than what is observed in the simulations. Despite 
these slight discrepancies, which are usually about 5% and no larger than about 10%, we 
consider this remarkable agreement given the simple monochromatic trial function and the 
linearization in the perturbation treatment. Figure 0] shows the details of the perturbation 
variables as calculated from the theory, and they reveal that the actual oscillation frequency 
uj = uo s + 58 is only insignificantly different from uo s . Direct comparisons between the 
predicted modulation Q and the simulated values as a function of microwave amplitude for 
different values of rj and a are summarized in Figures E] and El Comparisons for Qr are 
provided in Figure El and the attenuation f3 is validated in Figure El The general agreement 
is very good for all parameters tested. We observe the simple relationship (3 = a/2 for 
all oscillating solutions and the general trend of the polynomial (Equation (fTH|) or p7jl ) 
is to provide a modulation frequency slightly larger than what is observed in simulations. 
This is quantitatively similar to the values obtained from the energy balance perturbation 



theory developed in Ref. |17|, where the true modulation frequency was consistently slightly 
under-estimated. We finally observe the broad applicability of the asymptotic value of the 
modulation frequency given in Equation f)27j) . 

V. RABI-TYPE OSCILLATIONS 

As was argued in Ref. [17(, the Rabi-oscillations observed in, e.g., Refs. may 
be closely related to the classical transient modulations described above. The numerical 
simulations in Ref. [l^ were conducted in close agreement with the procedures described in 
Ref. I2I, by applying a small temperature T « 50mK (kT/Hj ~ 10~ 2 C 1) to a Josephson 
junction, which is perturbed by a microwave pulse starting at t = with a frequency 



lu s = \J\ — r\ l . Varying either the microwave amplitude or duration may yield Rabi-type 
oscillations in switching probability when a subsequent microwave pulse with a smaller 
frequency u p < uj s is applied to probe the energy state of the junction as it is left by the 
first microwave pulse. The frequency of the second (probe) field was chosen to match an 
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FIG. 3: Normalized time response in phase amplitude A and average phase ipo (a and c) and 
normalized energy H — Hq (b and d) as a function of time after onset of microwave perturbation. 
Figures a and b show results of the perturbation theory Equations (f7jl-pij) and ()28[). Figures c and d 
show results of direct simulations of Equation (^Q) for the same parameters as (a and b). Parameters 
are: a = 0.00151477, r) = y/l - u% = 0.904706, e s (t) = 0.00108 • 0(t), and T = e p (t) = 0. 

anharmonic amplitude that could excite the junction beyond the energy saddle point and 
lead to escape from the well. The measurements were conducted for varying microwave 
amplitude and fixed microwave duration. 

The experiments described in Ref. |5( were conducted slightly differently. First, the 
Josephson system was a small-inductance (3l interferometer, which we can conveniently 
approximate with a single Josephson junction, since the dynamics for small /3l is known 

|. Second, the probe pulse was 



to be well represented by a single degree of freedom 
in this case a "dc" pulse as sketched in Figure ^ and the measurements were presented 
for fixed microwave amplitude and variable duration. We will here provide simulations of 
the classical system described in Equation ((TJ and parameterized by information provided 



in Ref. 0. Characteristic current and frequency are I c = 6.056fxA and cu « 110 • 10 s 
which lead to a normalized temperature of kT/Hj w 2 • 10~ 4 (T = 30mK), and a normalized 
microwave frequency u s = 0.6527147, which is close to the resonance such that rj — Jl — u^. 
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We have estimated the dissipation parameter a from the reported decay of Rabi-oscillations, 
a = 0.00151477. The normalized applied microwave amplitude e s is varied between and 
0.01, and the normalized duration is in the range [0; 3000]. Simulations are conducted with 
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FIG. 4: Linear transient responses to microwave onset versus normalized time as calculated from 
Equations ()28 [) -(|3U [) for the parameter values of Figure 01 

a probe pulse s p (t) as shown in Figure followed by a short time in which we determine 
whether or not the junction has switched from the zero-voltage state. For every specific 
set of parameters, this type of simulation is conducted 25,000 times, each with a randomly 
chosen value of microwave phase 9 S , and the switching probability is recorded, before a new 
microwave duration is chosen. Typical results are displayed in Figured, where the switching 
probability P, reported as population occupancy of excited quantum state in Ref. 0, is 
shown as a function of microwave duration. We clearly observe the Rabi-type oscillations of 
this classical system with a distinct frequency. Moreover, the oscillation frequency is in near 
perfect agreement with the reported comparable figure of the experimental measurements 
as well as the theoretical value Qr provided by the analysis in the above section. Figure UJp 
shows the simulated system energy as an average of many thermal realizations, each with a 
randomly chosen 9 S , of the trajectory (at T »s 30mK). The different choices of 8 S G [— tt; it] 
yield slightly different trajectories of 6A(t), S6(t), and 8(po(t), which in turn make possible 
switching a function of 9 S . The importance of this ensemble average can be viewed in 
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FIG. 5: Normalized modulation frequency £Ir (at T = 0) as a function of normalized signal 



amplitude e s for (a) a range of damping a values, with r) = 0.94259 (w s = yl — rj 2 = 0.577886). 
Gray line with large "+" markers indicates the asymptotic form given in Equation (J27|). Arrows 
indicate the frequency for which small oscillations would be overdamped for respective values of 
a. (b) A range of dc bias rj values, with a = 0.001. Curves represent solutions to Equation (|17|) . 
Markers represent data from numerical simulations of Equation Q with e p (t) = T = 0. 

Figure Et, where a single trajectory of energy is shown for the applied temperature of T = 
30mK. After a short transient, the thermal effects are dominating the individual trajectory, 
exciting the modulation frequency at random phase. Only after an appropriate average 
over trajectories do we observe the smoothly decaying envelope of the modulated energy 
curve in Figure [7)d, which in turn is very similar to the (single) zero-temperature trajectory 
seen in Figure [7|i. The close connection between the transient energy modulation and the 
Rabi-type modulation is obvious, and the close relationship between classical numerical 
simulations, classical theory, and the experimental measurements is further emphasized in 
Figure |H1 where the experimental data (open markers) of Ref. || are shown alongside theory 
(solid curve) and simulations (closed markers). The three-way close agreement indicates 
that much about these measurements can be understood from classical theory of driven, 
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FIG. 6: Normalized attenuation (5 as a function of damping a and normalized signal amplitude 
e s for T = 0. Lines represent the derived relationship (3 = a/2 and markers represent data from 
the simulations of Figure 03 

damped Josephson junctions. 

VI. CONCLUSION 

We have presented a perturbation analysis of the classical, nonlinear model describing the 
microwave-perturbed Josephson junction. Results show direct quantitative analogy between 
experimentally reported Rabi-oscillations and the classical transients to phase- locking, and 
results are presented as specific functions of experimental system parameters. The analysis in 
this paper is both a complement and an extension to the analysis of Ref. where an energy 
balance perturbation approach was applied to produce similar connection between reports of 
Rabi-oscillations and the nonlinear classical Josephson system. The extension provides direct 
closed form solutions of a (Rabi-type) modulation frequency as well as attenuation, and the 
resulting frequency agrees closely with the previously obtained result. This consistency lends 
credibility to the value of a classical interpretation of the reported Rabi-oscillations, since 
the quantum mechanical and classical intuition turns out to be very similar (as is also the 
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FIG. 7: (a) Escape probability P as a function of signal duration AtgUjQ 1 . Numerically simulated 
quantities are represented by dots (each of which represents 25,000 escape events), (b) Ensemble 
average of normalized energy, simulated through Equation as a function of signal duration 
(ensemble size N = 50,000, randomly chosen 8 S ). Parameters are: T = 30 mK, a = 0.00151477, 



T] = y/l-v* = 0.904706, lo s = 2itv 01 /uj = 0.652714, e s = 0.00217000, e p = 0.0847400, and 
loo ~ 110- 10 9 s _1 . These parameters are inspired by those reported in Ref. [5j. (c) Single trajectory 
of normalized energy versus time for parameters listed for (b). (d) Single trajectory of normalized 
energy versus time for T = 0. 



case in laser physics 24]). 

The macroscopic quantum picture interprets the observed oscillations as a result of a mi- 
crowave induced temporal variation in the population probability of two, or more, intrinsic 
quantum mechanical energy levels in the Josephson washboard potential. This slow vari- 
ation is then indirectly observed through the associated variation in tunneling probability, 
switching from the zero-voltage state as a result of the application of the probe pulse. The 
stochastic nature of the system is due to the absorption of microwave photons as well as the 
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FIG. 8: Modulation frequency f2/j as a function of normalized signal amplitude e s for a = 
0.00151477, 77 = 0.904706, T = 30 mK. Lines represent calculations using Equation (|17|) . o represent 
statistical data from simulations, and • are the measurements copied from Ref. |^|. 

subsequent tunneling probability of the measurement. 

The classical picture presents a system with a microwave induced temporally modulated 
energy. The slow variation is indirectly observed through the associated variation in escape 
from the energy well when the probe is applied, with high probability for passing the en- 
ergetic saddle point during times of elevated energy and relatively little escape probability 
during times of depressed energy content. The stochastic nature of this system is due to the 
random phase of the microwave signal as well as the thermal fluctuations. 

Additional analogies between intrinsic quantum mechanical energy levels and the multi- 
valued resonances of the classical nonlinear system briefly outlined above (as well as in 
Ref. [17]) further bridges the gap between what can be expected from the two interpreta- 
tions of the microwave induced measurements. Since the experimental measurements are 
concerned only with detecting the escape event, and since the classical and quantum mechan- 
ical interpretations seem to provide the same signatures for that escape, we are faced with a 
fundamental ambiguity of how to read this information. Potential applications of Josephson 
technology for quantum information processing will therefore benefit from the development 
of new unambiguous measurements, which must present signatures of macroscopic quantum 
behavior that cannot be explained classically. However, we submit that, given the very close 
quantitative agreement between our classical theory and the experimental measurements, 
the intuition and the closed form expressions presented in this paper can be directly applied 
to guide future experiments and microwave manipulations of Josephson systems. 
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